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1. Introduction 

Traditionally the physical applications of integrable systems are exhausted by the 
approximations to real dynamical systems. The configurations of the (classical or quantum) 
model form the phase space which is a manifold M with the symplectic two-form uj or 
Poisson bi-vector tt. 

The phase space might carry a natural complex structure, such that the symplectic 
form w is a holomorphic (2, 0)-form, the Hamiltonians H{p, q) are the holomorphic func- 
tions and the vector fields are holomorphic vector fields (see [|l| for example). The use of 
integrable systems as describing the evolution in the physical models is less transparent in 
this case. 

The integrable systems in the holomorphic sense entered physics approximately at the 
same time as the string theory did. Particle which has a one-real-dimensional worldline is 
naturally described with the help of a real phase space, its (real) time evolution being gen- 
erated by the real Hamiltonian. The worldsheet of a string is a complex curve, embedded 
into the target space. It senses holomorphic geometry in many different ways. In partic- 
ular, the embeddings of the worldsheet E are governed by a two dimensional conformal 
field theory on E. An example of holomorphic integrable system relevant to the latter is 
the famous Hitchin system. The phase space of this model is the cotangent bundle to the 
moduli space 7V1g(E) of holomorphic G-bundles over a Riemann surface E. One can think 
of Knizhnik-Zamolodchikov-Bernard [|| equations in two dimensional WZW theories as of 
the non-stationary quantum version of Hitchin system @] IQ] 10] • The role of times 
is played by the complex (and hypothetically W)-moduli of E. Complex time evolution 
occurs also in the models of A^ = (2, 2) strings, where space-time may have (2, 2) signature 

i- 

However there exist other possibilities for integrable system to encode the physical 
information. 

In particular, the rich source of holomorphic integrable systems is the combination of 
supersymmetry and duality. It is known for some time now that the holomorphy of certain 
quantities (like the superpotential in A/" = 1 or prepotential in A/" = 2) in the supersym- 
metric theories in three/four dimensions yields powerful predictions for the behavior of the 



quantum theory even in the presence of the non-perturbative effects |jT^, |jTT[], [p!2| , [p!3| , 
P^ , [|T5| . In particular, the complex structure of the moduli space of vacua in A/" = 4 
3(i gauge theories and N = 2 Ad gauge theories can be determined revealing the exciting 



link with the special nature of the geometry of the phase spaces of complex integrable 



systems [16], [O], |T3], [15|. The action variables appear as the central charges in the BPS 



representations of the susy algebra [16 



There exists an approach to a class of integrable system which allows to uncover the 
origin of their integrability/solvability. Namely, one realizes the system under investigation 
as a projection of a simple system on a larger phase space [jl^ , |2^] . This idea is actually 



a counterpart of the main principle behind the gauge theories - the complex dynamics 
of the actual world (as far as most of the fundamental interactions are concerned) is a 
projection of a somewhat simpler dynamics of the extended phase space. One of the goals 
of the present discussion is to use the analogy between the two ideas and explain certain 
properties of integrable systems as well as gauge theories. 

We are going to study the phenomenon of duality whose precise definition is presented 
shortly. Duality is a subject of much recent investigation in the context of (supersymmetric) 
gauge theories, in which case the duality is an involution, which maps the observables of 
one theory to those of another. The duality is powerful when the coupling constant in one 
theory is inverse of that in another (or more generally, when small coupling is mapped 
to the strong one). For example, a weakly coupled (magnetic) theory can be dual to the 
strongly coupled (electric) theory thus making possible to understand the strong coupling 
behavior of the latter. In particular, it was shown by N. Seiberg and E. Witten [^ that 



using the concept of duality one can find exact low-energy Lagrangian of A/" = 2,(i = 4 
SU{2) gauge theory. A more fascinating recent development is that the duality connecting 
weak and strong coupling regimes of one or different theories may have a geometric origin. 
The most notorious example of that is provided by M-theory [|2l|,[|2^. We are going to 



study the dualities in integrable systems, related to the gauge theories with the emphasis 
on their geometric origin. 

The study of geometry of integrable systems also allows to understand the origin of 
certain constructions of separation of variables [^ . The similarity of this construction to 



the description of the D2-brane moduli space and its role in the understanding the string 
duality makes one hope that both subjects - many-body integrable systems and gauge 
theories (more generally, D-geometry of M. Douglas [^[^) will benefit more from each 
other in the near future. 

Topics left beyond the scope of the paper. To keep the size of the paper within 
reasonable limits we decided to restrict our attention with the pure many-body systems. 
More or less everything we have said can be carried over to the spin systems both of the 



'adjoint' |2^|^ and 'fundamental' [^ type. We don't discuss extensively the relation 



of our dualities in integrable systems to the physics of D-branes [^[0[|3T|]. Some of 



the results in this direction together with the applications to the theory of separation 
of variables can be found in ||2^. Also, except for the general discussion and two-body 



examples we don't treat quantum case. For some results related to our main topic see 
|3^|^^[Q[^. Realizations of elliptic Ruijsenaars-Schneider models via Hamiltonian and 



Poisson reductions can be found in ||36|| . 



Organization of the paper. Various concepts of duality are discussed in the section 2. 
The examples of the dual systems are studied in section 3 where mostly two-body case is 
treated, both classical and quantum one. Many-body systems are studied in the section 
4 with the explanation of the dualities between them coming from Hamiltonian/Poisson 
reductions. The section 5 is devoted to the gauge dynamics and their relation to the 
integrable systems discussed so far. We discuss the geometry of the moduli spaces of vacua 
of supersymmetric gauge theories in three, four, five and six dimensions and construct a 
little dictionary translating the notions of integrable systems to those of gauge theories. 
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2. The concepts of Duality: 

Let (M, u) be a symplectic manifold. There exist Darboiix local coordinates (cf. 
in which the symplectic form looks like a canonical one: 

m 

iv = ^dpiAdq^ (2.1) 

The local canonical coordinates are defined up to the symplectomorphisms. Unlike the gen- 
eral diffeomorphisms, which have N functional degrees of freedom, N being the dimension 
of the manifold, the symplectomorphisms have only 1 functional degree of freedom. 

The evolution of a Hamiltonian system is defined with the help of Hamiltonian H : 
M — i> IR. The function H defines a Hamiltonian vector field by the formula 

iy^u = -dH (2.2) 

The integrable system on M has a maximal collection of the functionally independent 
commuting Hamiltonians if^ z = 1, . . . , m = |dimM : 

[Vh.,Vh,] = Q (2.3) 

Let h : M — > S :^ IR"^ be the map defined as: h : a; i-^ {Hi{x), . . . , Hm{x)). Liouville's 
theorem states that the integrable system has a normal form locally: there are coordinates 
(/i, (/?*), such that 

u = y2 dli A dip'' 

(2.4) 

Hk = fk{{I}) 

i.e. Ik are coordinates on B. For a sufficiently small domain U G B the space h~^{U) is 
the product f/xR"""" xT^"'"'' and (p^ are standard linear coordinates on IR''""' xT^"""". 
If the common level set of all Hamiltonians is compact then this set is isomorphic to the 
torus of the dimension m. In that case one may impose a condition on the coordinates 
ip^ that the differentials d<p>^ have periods which are integer multiples of 27r. This fixes the 
coordinates (/, ip) up to the action of discrete group PGLm(2). 

The Liouville theory also has a counterpart in the holomorphic setting where the 
manifold M is replaced by the complex manifold, the symplectic form is a holomorphic 
closed (2, 0)-form, the Hamiltonians H are the holomorphic functions and the vector fields 
are holomorphic vector fields. The Liouville theorem modifies in this case. In fact the 
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Liouville real tori are replaced by the complex tori. If we require these tori to be abelian 



varieties then we get what is called algebraically integrable system [^. In the family of 
such varieties the degenerate fibers can appear. 

The coordinates li are referred to as "action" variables. If n = m then there is a nice 
formula for /. Let bi,b2 E U G B he sufficiently close to each other. Choose a basis e^ in 
IHi(h~^(6i), Z). Connect the points 61 and 62 with a path 7 C t/. The base 6^^ can be 
transported to lHi(h~^(7), Z) by means of the Gaufi-Manin connection and it defines an 
element f G H2(h-i(7), h'^ibi U 62); 2). Then 

I{b2) - I{bi) = fuj 

2.1. {p,q)-^{I,ip) 

Suppose that we have two integrable systems {Hk} and {H^} on the same symplectic 
manifold M . In this situation we say that these two systems are dual to each other. 

Notice that this definition does not make duality an involution. 

A pair of integrable systems given on one symplectic manifold (M, u) is called self-dual 
if there exists a symplectic involution a : M ^ M exchanging {Hk} and {H^}, i.e., such 
that for any k — 1, . . . ,n 

a*Hk = H^ 



Once we have two integrable systems on the same manifold such that both collections 
of Hamiltonians constitute at least locally a coordinate system on the phase space M, one 
can write down the equations of motion of the second integrable system in the second 
order formalism using the action variables li of the first system as the coordinates q'^ for 
the second. 

The global version of this definition is: Two Hamiltonian systems are dual to each 
other in the sense of action-coordinate (AC) duality if the action variables li of the first 
system coincide with the coordinates g* of the second one and vice versa. 

2.2. I ^I^ 

In the holomorphic algebraic category there is an interesting complication: the torus 
has a complex dimension m and therefore lHi(Tb; Z) = 2 "^. 

Again, let 61, 62 G t/ C -B be sufficiently close to each other. Choose a symplectic 
basis Cb, = (Aa,S^) in Hi(h-i(6i), Z) such that A n A = B n B = 0, A^ n B^ = 6^, 



where fl is an intersection form Hi ®-^ Hi -^ Z. Connect the points hi and 62 with a path 
7 C t/. The base eb^ can be transported to Hi(h~^(7), Z) by means of the Gaufi-Manin 
connection and it defines an element F G H2(h~^(7), h~^(6iU62); ^)- The (A, i?)-sphtting 
of the basis e^^ is also transported via Gaufi-Manin along 7 and defines the decomposition 
f = Ta + f B, f A,B e H2(h-i(7), h-i(6i U 62); 2). Then define 



1(62) - /(61) = / 

P{h2)-P{hi)= [ U 
Jfn 



A 



(2.5) 



We can take as the action variables the components of /. The reason for the S-cycles to be 
discarded is simply the fact that the S-periods I^ are not independent of the A-periods. 
On the other hand, one can choose as the independent periods the integrals of u over any 
Lagrangian (in the sense of fl) sublattice in Hi(Tb; Z) transported along 7. 

This leads to the following structure of the action variables in the holomorphic setting. 
Locally over a patch in B one chooses a basis in Hi of the fiber together with the set of 
A-cycles. This choice may differ over another patch. Over the intersection of these discs 
one has a Sp2^(Z) transformation relating the bases. Altogether they form an Sp2^(Z) 
bundle. It is an easy excercise on the properties of the period matrix of abelian varieties 
that the two form: 

dP A rf/f (2.6) 

vanishes. Therefore one can always locally find a function JF - prepotential, such that: 

This duality maps the integrable system to itself. It is called action-action (AA) dual- 
ity. 

2.3. Quantum, duality 

There exists a clear quantum counterpart of this picture. Consider the eigenvalue 
problem for the Schrodinger operators and the issue of the normalization of the wave- 
functions. 

The quantum integrable system is a complete collection of "independent" (in the 
appropriate sense) commuting operators {Hi},i = l,...,rM, acting in the Hilbert space 



7i of the model. By completeness we mean that these operators have simple common 
spectrum U C IR"^.The commuting operators have common eigenf unctions. Generically 
the eigen- value problem: 

Hk\X) = ek{X)\X) (2.8) 

has the unique (up to normalization) solution. Here e^ is the corresponding eigenvalue 
and A is a label, which takes values in some set A. Altogether, Cfc form an imbedding 
e* : A ^ [/ 

Typically one has another set of commuting operators ("position operators") {H^} 
and the eigen-states \x) with eigen-values efc(x) in the Hilbert space of the model are 
represented as the appropriate functionals on the space of the eigen-values of the operators 
Hj^. Here x is another label, which takes values in the set A''^, which is mapped by ef 
to U^ e IR"". The famihar case is M = T*M, Hm = L'^{M), the operators Hi are 
represented as commuting differential operators, and H^ are represented as operators of 
multiplication by a function ef{x). Suppose that we are given two classically AC dual 
Hamiltonian systems. Let li and //^ be their action variables (recall that //^ are the 
coordinates for the first system). Assume that there exist quantum integrals of motion for 
both systems. Let us denote by li the quantum integrals of motion of the first system and 
by I^ the quantum integrals of motion of the second one. 

Once we have a quantum integrable system we can identify Ti with the space -^^(A) 
of square integrable functions on A (w.r.t. the spectral measure dn). Indeed, choose a 
basis in 7i consisting of the common eigenvectors ( p.8| ) |A), where A = {Ai, . . . , A^} G A. 
Then for any |'i/') G Ti one can associate the function {\\ip) on A. Of course this mapping 
Ti -^ -^^(A) depends on the normalization of the eigenvectors we have chosen. 

In particular any operator A acting on Ti can be expressed as an operator acting on 
L2(A) as 

i:V(A)^ ! (X\A\\')i:(X')d^i{\') 



JK 

Now suppose that we have two integrable systems Hi , . . . , H^ and H^ , ■ • • , H^ on the 
same Hilbert space Ti. We can use the first one to identify the Hilbert space with the 
space of functions on its spectrum -^^(A) and write down the operators of the second 
integrable system as acting on these functions and not on some abstract Hilbert space 
vectors. Consider the function (A|A^) G -^^(A) ® L'^{A^) = L'^{A x A-^), where |A) and 
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I A^ ) are the eigenvectors of the first and the second integrable system respectively. This 
function by definition satisfies for any k — 1, . . . ,m the foUowing equations: 



(Al^f |A')(A'|A^)d^(A') = ef (A^)(A|A^) 

A 

(A|A^')(A'"|^.|A^)rf;u(A^') = e,(A)(A|A^) 

AD 



(2.9) 



We see that the function (A|A^) turns out to be an eigenfunction for two commuting set of 
operators acting on two groups of variables. Otherwise the function (A|A^) is not uniquely 
defined by the two dual integrable systems since the arbitrary change of the normalizations 
of the eigenvectors |A) i-^ F(A)|A) and |A^) t-^ F^ {X^)\X^) results in 



(A|A-^)^F(A)F^(A^)(A|A^) 



Note also that though we have written the equations (|2.9[ ) as integral equation with 
the kernel being a generalized function for particular dual systems these equations may be 
differential or difference ones. (In fact it happens in examples considered in the sequel.) 

Analogously to the classical case a pair of quantum integrable systems {Hi} and 
{Hf} is called self-dual if there exists a unitary involution a : Ti ^ Ti exchanging the two 
collections of operators, i.e. such that for any /c = 1, . . . , m 

a-Hk=H^-a 

Though the two integrable systems of a self-dual pair are completely equivalent, the corre- 
sponding function (A|A'^) does not necessarily satisfy the condition (A|A^) = (A^|A). One 
can make it obey this equation after a suitable normalization of the eigenvector bases. 

In some cases it is natural to choose A among the action variables of the classical 
integrable system. 

3. Examples of dual systems. One degree of freedom: 

In this section we work out explicitly a few examples of the dual systems. 
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3.1. Classical systems 

Two-particle systems which we are going to consider reduce (after exclusion of the 
center of mass motion) to a one-dimensional problem. The action-angle variables can be 
written explicitly and the dual system emerges immediately once the natural Hamiltonians 
are chosen. The problem is the following. Suppose the phase space is coordinatized by 
(p, q). The dual Hamiltonian (in the sense of AC duality) is a function of q expressed in 
terms of /, </?, where /, </? are the action-angle variables of the original system : H£i{I, tp) = 
Ho (q) . In all the cases below there is a natural choice of H^ (q) . 
Calogero oscillator. The Hamiltonian in the center of mass frame reads as: 



2 2 2 2 



(3.1) 



where u) and v are the parameters. In the limiting cases u = and a; = one gets the usual 
oscillator and rational Calogero-Moser system respectively. The action-angle variables /, (f 
can be found by the standard procedure: 



with the result: 



hf'"''' = i 



'2(^-^-$M. 



dq(dl\-^ 



E — uju 1 



Hd{I.^) 



2u Au 

q^ I 



p" + {uq - -f 

q 



l + - + Jl + yCosv. 



(3.2) 



(3.3) 



2 u 

The limit z/ -^ is straightforward, yet tricky. We must rescale </;—!► 2 </? since the period 
of motion jumps as z/ approaches zero. We get: 



E 
2^ 



Hd{I, 'f)=q = 2VIcos{lp) 



(3.4) 



The limit a; ^ is more subtle as the classical motion becomes infinite. For the system 
with the Hamiltonian 

HiP^^) = Y + ^ (3.5) 
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the action variable could be defined as the asymptotic value of the momentum: / = \f2E. 
This choice gives rise to the evolution, linear in the "angle" -like variable, 



Z/2 



2 2 2 

Vf IT \ ^ '^\^ 



Sutherland model. The Hamiltonian is: 



1 , 



The action variable / can be chosen to be: 

I^V2E (3.8) 

To prove that one might go to the coordinate t = cos(q') and compute the integral -^ § pdq 
by residues. The angle variable 99 can be determined from the condition dpAdq = dl Ad<f>. 
We get: 

d,= ^i^ (3.9) 

sin*^(g) 



2iy2 
i^D (/,(/?)= cos(g) = cosv?Wl -J- (3.10) 

Notice, that ( p.lO|) coincides with the Hamiltonian of the rational Ruijsenaars model (see 

below) . 

Elliptic Calogero — Moser system. The Hamiltonian is: 

2 
H{p,q) = ^ + u^Pr{q) (3.11) 

Here p, q are complex, pr(<?) is the Weierstrass function on the elliptic curve Et-: 

-^ 1 1 

Priq) = — + Y] ^- : 7^-7 : r^ (3.12) 

(to, n) G "^ 
(TO,n)^(0,0) 

Let us introduce the Weierstrass notations: x — priQ)^ y = pr('?) ■ We have an equation 
defining the curve Er'. 

3 3 

y2 = 4x3 - g^(r)x - gs{r) = 4 J](x - e,), J^ e, = (3.13) 
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The holomorphic differential dq on E^ equals dq = dx/y. Introduce the variable cq = E/u'^. 
The action variable is one of the periods of the differential ^;^ on the curve E = H(p, q) : 



I=^i dq^ 2{E-v^pM) = ^ f , ^"'^^^^^ = (3.14) 

271- J A 2 V27r J A V (x -ei){x- 62) {x - 63) 

The angle variable can be determined from the condition dp A dq = dl A dip: 

d^= c. -rrt F^ / = (3-15) 

where T{E) normalizes dip in such a way that the A period of dip is equal to 27r: 

T(E) = -L / '''• (3,16) 

Thus: 



2iT{E)dp) 



udip 



4^^=o(^-ei 



4nli(^-i.: 



where 






-2iT{E)y/ 601602603 =— f 



+ -E- *. = -E 



a; - eo 3 ^ eoi 3 ^ eojCoj 

1=1 7=1 



C^j 02 Oj 



Introduce a meromorphic function on Et-: 



(3.17) 



1 1^1 Iv^ e,, (3.18) 



^Jz) = . r--^ (3.19) 

V x-es 

where z has periods 27r and 27rr. It is an elliptic analogue of the cosine (in fact, up to a 
rescaling of z it coincides with the Jacobi elliptic cosine) . Then we have: 



HD{I,^)=cnr{z)=cfir„{^)\\l- —=, — ^A^— (3.20) 

y ztj — z/^ea 
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where te is the modular parameter of the relevant spectral curve f^ = 4]^j^j^(t — ti): 

dt \ // /■ dt 



Te 



B 



m=iit-u 



A 



mt^s-u 



For large /, 2E{I) ~ P 

Elliptic Ruijsenaars model. The Hamiltonian is: 

H{p, q) = cos(/3p)v/l-2(/3z.)2p,(g). 



(3.21) 



(3.22) 



As the curve Ej. degenerates one flows down to the trigonometric (Pt(?) -^ ^.^L \ ) or 
rational (Pt(q') -^ A-) Ruijsenaars system. The spectral curve H{p, q) = E helps to define 



the action variable /: 



up to the transformations 



27r J A 



(3.23) 



rD 



27r 



where ni,n2,n3 G 2 and (ni,?7,2) = 1. These transformations reflect the freedom in the 



choice of the cycle A on the spectral curve. The appearence of ns was used in [^ . We 
can write an explicit formula for the quantity which is better deflned: 



dl 



1 



dE 2V2.^v 7. /j^T^^^r; 



dx 



where now cq 



l-E-^ 



(3.24) 



. Under A -^ B transformation 4^ gets multiplied by te, where te 



is deflned as in ( p.21| ) . Quite similarly to (|3.17|) we get: 



dip 



1 



dx 



^(^'x/fiu^^ 



(3.25) 



with 



T{E) = — 



27r 



dx 



A 



U^=oi^~ei) 



(3.26) 



Finally, for Hd given by (|3.19| ) we get: 



,. — /^ — /^/i 2(/3z^) ei3 

Hd U,f) = CUr (Z) = CUrp. ( V^) W 1 9 



(3.27) 



14 



Asymptotically, for large /, E{I) ~ cos{(3I). 

General elliptic model. In the general case one modifies the formula ( |3.13| ) in such a 
way that the coefficients g2 and g^ are the sections of the line bundles 0{4n) and 0{6n) 
respectively over B ~ P . The elliptic curve E^^ defined by the modified ( |3.13[ ) degenerates 
over the divisor of zeroes of its discriminant: 

A = gl- 27gl (3.28) 

which is a section of 0{12n). The latter has generically 12n zeroes. To make the total 
space of fibration isomorphic to the K3 surface (compact simply-connected symplectic 
surface) we need 24 singular fibers, i.e. n = 2. The Hamiltonian of the integrable system 
we consider is any function on B. It gives rise to a meromorphic vector field on M, which 
linearizes along the elliptic fibers. The symplectic form is given by: 

dx Adz , , 

LU = 3.29 

y 

where z is the projective coordinate onB. Under change of the variables: z = -, y = — -^, 
X = ^ the form ( p.29| ) goes over to ^^4^ and the equation ( |3.13| ) is mapped to 

y2 = 4x^ - g2(z)x - gsiz) 

where the polynomials gk are defined through the relation: 

gkiz)=z'"'gk{l/z). 

Over a simply connected region C7 C P \A~^(0) one can trivialize the bundle of the 
first homologies lHi(i?2, 2), in particular to make a well-defined choice of the A-cycle of 
the elliptic fibre. The local action variable / = I{z) is defined over U by the equation: 

^ = T(z) = ^ / ^ (3.30) 

dz 2n J A y 

where the integral is taken over a chosen A-cycle. The fibration of lHi(i?;j,Z) over 
B\A~^(0) is non-trivial and there is no global monodromy invariant choice of A-cycles. 
So the action variable is defined by ( |3.30|) only locally. The monodromies around the de- 
generate fibers corresponding to various singularities has been worked out by Kodaira and 



their physical interpretation can be found in [^. For generic polynomials (72(2;), (73(2;) the 
singularities are of the type Ai. 
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The angle variable dual to I{z) is nothing but the linear coordinate on the Jacobian of 
the fiber elliptic curve (|3.13|) . In particular it is periodic with the periods 27r and 2tvt{z). 



It is to be found from the relation: 

d^=^-. (3.31) 

T{z) y 

We can get a dual system by treating x as the Hamiltonian. Since x is not a mero- 
morphic function on K3 (it changes under the z ^ - transformation) this is only possible 
if we delete the elliptic fiber Eoo ■ 

Let us see what will be the action-angle variables. First of all, generically the fiber Cx 
over X eF is an incomplete hyperelliptic curve of genus 5. The holomorphic differentials 

on this curve are: 

z dz 
LOk = — F ^, k = 0,...4 

a/4x3 - g2{z)x - gs{z) 

The action variable I^ = I^{x) obeys the equation: 

"'^ r-(.) = ^/^ (3.32) 



dx 2% /^ y 

where L is a one-cycle in lHi(Ca;,2). Here we face AA duality in its extreme form: 
the freedom to choose L is much bigger then in the case of original system, since the 
corresponding duality group is Sp]^o(^)- We can partially integrate ( |3.32| ) to get 

/^ = — (/ ifT{z)dz (3.33) 

27r Jl 



where 



x = -— ^p(v?;r(2)) 



The angle variable is one of the linear coordinates on the Jacobian variety of C^, which is 

5 dimensional abelian variety: 

. D I dz 

The embedding of the Liouville tori into the abelian varieties of higher rank originating 
from hyperelliptic curves is a well-known phenomenon in the theory of integrable systems, 
going back to the original work of S. Novikov and A. Veselov [^ . 
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3.2. Quantum systems. 



Here we work out a few examples of quantum dual systems. 
Harmonic oscillator. The Hamiltonian (|3.1|) in the limit u = quantizes to: 



H 



1 d^ u^q^ 



2dr 



+ 



(3.34) 



Its normalized eigen- functions are |[42|| : 



Hi(jn =u;{n+ -)il}r. 






(3.35) 



HniqVuJ 



where Hn{$,) is the Hermite polynomial: Hn{^) = e^ (— 9^)"e~^ . Using this representa- 
tion of the wave-function one can easily obtain a reccurence relation (details are in the 
appendix) : 

Vn + It/jn+iiq) + ^/ni/^n-iix) = V2cJ#n(x) (3.36) 

It means that i(^n{q) is an eigen-function of the following difference operator: 

HD=T+y^+ v^T_ , T± = e± * (3.37) 



acting on the subscript n. It is easy to recognize in (|3.37|) the quantized version of (|3.4|) . 
Sutherland model. Here we deal with the Hamiltonian: 



H 



1 d^ v{v - 1) 

2'd^^ 2sin2(Q) 



(3.38) 



Its normalized eigen- functions are p2 

2 

n 






(3.39) 



For simplicity we take u and n to be half-integers. One can change u ^ —u — 1 to 
get another eigen-function with the same eigenvalue. Using the fact that the generating 
function for H?'s is 



z{y,x) = J2y'uf 



1 



1=0 



1 — 2xy + y^ 



(3.40) 
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one derives the recurrence relations (details are in the appendix): 

I + 1 — m „^ I + m. 

(3.41) 



^ 2/ + 1 '+^ 2/ + 1 ^"^ 



1 / l-z/(z/-l) , l-z/(z/-l) , 

cos^ V'n =-[\ J — ——tlJn+i + \ J T^^^-i 

2\y n[n + l) y n{n — 1) 

that is 'i/'n is an eigen-function of the finite-difference operator acting on the n subscript: 

HDi^iq) = cos{q)tjj{q) 

y n{n — 1) y n{n — 1) 

which is a quantum version of ( p.lO| ). 

Moral of the story. The moral of the previous discussion is that the polynomial de- 
pendence on momenta of the hamiltonian is traded with the rational potential of the 
dual system. The trigonometric potential is mapped to the trigonometric (= relativis- 
tic) dependence on momenta of dual system. The elliptic potential gives rise to elliptic 
(="double-relativistic" ) dependence on momentum of the dual system Hamiltonian. When 
the system with trigonometric dependence on momentum is quantized its Hamiltonian be- 
comes a finite- difference operator. The wave-functions become the functions of the discrete 
variables. The origin of this is in the Bohr-Sommerfeld quantization condition. Indeed, 
since the trigonometric dependence of momenta implies that the leaves of the polarization 
are compact and moreover non-simply connected the covariantly constant sections of the 
prequantization connection along the polarization fiber generically seases to exist. It is 
only for special "quantized" values of the action variables that the section exists. In the 
elliptic case the quantum dual Hamiltonian is going to be a difference operator of the infi- 
nite order. The self-dual elliptic many-body system is still to be constructed. It seems that 
to achieve this goal one needs a notion of the Heisenberg double for the central extension 
of the two dimensional current group |^ . 



Example of the prepotential. To illustrate the meaning of the AA duality we look at 



the two-body system, relevant for the SU{2) J\f = 2 supersymmetric gauge theory [|T6 

2 

H=^+K^ cos{q) (3.43) 

with A^ being a complex number - the coupling constant of a two-body problem and at the 
same time a dynamically generated scale of the gauge theory. The action variable is given 
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by one of the periods of the differential pdq. Let us introduce more notations: x = cos(g), 
y = ^^^^^ , u = ^. Then the spectral curve, associated to the system ( |3.43p which is also 
a level set of the Hamiltonian can be written as follows: 

y'^ = {x - u){x^ - I) (3.44) 



which is exactly Seiberg-Witten curve [|T2| as it was first observed in [^. The periods are: 



,D I I X-U 



(3.45) 

They obey Picard-Fuchs equation: 

\d^^ A{u^-1)) \ /^j=° 

which can be used to write down an asymptotic expansion of the action variable near 
u = oo oi u = ±1 as well as that of prepotential ( |2.7|) . The AA duality is manifested 
in the fact that near u = oo (which corresponds to the high energy scattering in the 
two-body problem and also a perturbative regime of SU{2) gauge theory) the appropriate 
action variable is / (it experiences a monodromy / -^ —I as u goes around oo), while near 
u = 1 (which corresponds to the dynamics of the two-body system near the top of the 
potential and to the strongly coupled SU{2) gauge theory) the appropriate variable is I^ 
(which corresponds to a weakly coupled magnetic U{1) gauge theory and is actually well 
defined near u = 1 point) []12|. The monodromy invariant combination of the periods ^A 



11^ -2T = u (3.46) 

(whose origin is in the periods of Calabi-Yau manifolds on the one hand and in the proper- 
ties of anomaly in theory on the other) can be chosen as a global coordinate on the space 
of integrals of motion B. At w ^ oo the prepotential has an expansion of the form: 



T ~ -u log w -i- . . . ~ /^ log / -i- y^ — / 



2-4n 

n 
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3.3. Appendix. 

To derive the recurrence relation for the oscillator wave-functions we use the creation 
operator representation: ip^ = ~^(~^C + OV'n- Applying this relation twice and using 
the fact that i(^n is an eigen- function of H one arrives at ( p.36| ) . For the Sutherland model 
we use two obvious relations: 

ix-y)d,Z = ydyZ (3.47) 

(1 - 2xy + y^)dyZ = (x - y)Z (3.48) 

Next, ( |3:i7D implies: 

(ydy - m)ar^ = (x - y)d:^+'Z (3.49) 

and ( F48| ) yields: 

((1 - 2xy + y^)dy + y- x)d^Z = m(l + 2ydy)d^-^Z (3.50) 

Combination of those two gives rise to ( p.41|) . 
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4. Duality in Many-Body Systems: 

In the previous sections we discussed the concepts of duality and worked out explicitly 
several examples of dual two-body systems in both classical and quantum cases. We now 
turn to a study of many-body systems. The many-body systems can be divided into 
three classes: rational, trigonometric and elliptic ones. The Hamiltonian of the model 
may depend on momenta/coordinates in any one of these three fashions. The duality 
transformation exchanges them. 

4.1. Examples. 

We summarize the systems and their duals in the following table: 

rat. CM ^ rat. CM 
trig.CM ^ rat.RS (4.1) 

p^o T T R^Q 

trig.RS ^ trig.RS 



Here CM denotes Calogero-Moser models |^|^[0 and RS stands for Ruijsenaars- 



Schneider [Q [^ [|0l |J5T|] |]52| . The parameters R and /3 here are the radius of the circle 
the coordinates of the particles take values in and the inverse speed of light respectively. 
The horizontal arrows in this table are the dualities, relating the systems on the both 
sides. Most of them were discussed by Simon Ruijsenaars |J5^,|^. We notice that the 



duality transformations form a group which in the case of self-dual systems listed here 
contains SL2(Z). The generator S is the gorizontal arrow described below, while the T 
generator is in fact a certain finite time evolution of the original system (which is always 
a symplectomorphism, which maps the integrable system to the dual one) . We begin with 
recalling the Hamiltonians of these systems. Throughout this section q^j denotes qi — qj . 
Rational CM model. The phase space is {T*V)/r, where V is a linear space acted on 
by a Coxeter group F. We consider the simplest case V = M, , F = Sn+i. Let {pi,qi) 
be the set of coordinates, i = 1, . . . , A^ + 1 with the constraint J^Qi — J2Pi — ^- The 
Hamiltonians can be conveniently packaged using the Lax operator: 



1 

ffc 



m = ^TiL^ 



^ ■ (^ R \ (4-2) 

Lij =Pidij H '— 

Qi~ Qj 
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In particular, the quadratic Hamiltonian reads: 

,2 



Trigonometric CM = Sutherland model. The phase space is {T*V)/r, where V is a lin- 



7V-1 



ear space acted on by an affine Coxeter group f . We consider the simplest case V = 1R 
r = Sn X ^Z . Let {pi, qi) be the set of coordinates, i = 1, . . . ,N with the constraint 
^ q'i = ^Pi = 0, and the identifications Qi ^ Qi + ^ni.rii G Z. The Hamiltonians can be 
conveniently packaged using the Lax operator: 

^ ^^z^a-5^ (4-4) 

'' "^^^^ 2sin(fc^) 

In particular, the quadratic Hamiltonian equals: 



i<j 



4sin2(^i^yM) 



Rational RS = Relativistic rational CM model . The phase space is {T*V)/r, where 
y is a linear space acted on by an affine Coxeter group F. We consider the simplest case 
V = IR ~ , r = Sn X ^Z . Let {pi, qi) be the set of coordinates, i = 1, . . . , A^ with 
the constraint ^Qi = "^Pi = 0, and the identifications Pi ^ Pi + ^ni, rii G Z. The Ha 
miltonians can be conveniently packaged using the Lax operator: 

1, 



k 



/S" TT , (/S") 



'^''-^-""^^,Xl\^- 



2 (4.6) 



In particular, the Hamiltonian \{Hi — H^i) equals: 



if = llY(L + L-i) = 5^cos(/3p,)nJl-^- (4.7) 



J^^ V % 



The Lax operator ([4.6| ) is gauge equivalent to the operator 



^'^ '' q., + /3Z. V ^^ ^J 



-^PPiUt. u+, 



^^ ^ ^ ^uw{q^) - n(^) = n(^-^- 



(4.^ 
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In the limit (3^0 both L, L of (|]|),(U) behave as Id-z/? ( Lax operator in (|3) ) +o(/3). 
Trigonometric RS = Relativistic Sutherland model. The phase space is (T*y)/r^, 
where V is a linear space acted on by a double afRne Coxeter group Te, E being an 
elhptic curve. We consider the simplest case V = 1R^"\ T = Sn y< (^2^ © ^2^). Let 
{pi, Qi) be the set of coordinates, i = 1, . . . , A^ with the constraint J2^i — "l^Pi — ^^ ^^d 
the identifications Pi ^ Pi + ^^i, Qi ^ Qi + ^"^i^ ni^rrii E 2. The Hamiltonians 

can be conveniently packaged using the Lax operator: 



Hk 



k 



RPu 



^IJ 



M^iQ^j+H)t}j\ 



sm 



RPu 



sm 



Rq 



jfe 



In particular, the Hamiltonian h{Hi — -ff_i) equals: 



(4.9) 






s\n^{R[5v) 



sm 



2lR(q 



l iiK<li3) \ 



The Lax operator ([4. 61) is gauge equivalent to the operator 



sinf^ 






^± NR P(q,±Pu) ^, , -/t fR^ 



V ^ 



(4.10) 



(4.11) 



In the limit i? ^ 0, with (3 fixed the expressions ( [4.9| ),( ^T0|) ,( [4.11|) naturally go over to 
(|]|), (F^), (|48D respectively. In the limit /? ^ 0, i? fixed both L, C behave as ld-i(3{ Lax 
operator in ( ^ ) ) + o(/3). 



^.^. Explanations: Hamiltonian and/or Poisson reduction 

Suppose we are given a symplectic manifold {X,ujx) with the Hamiltonian action of 
a Lie group G with equivariant moment map p : X —>■ g* . The symplectic quotient of X 
with respect to G is the symplectic manifold M, denoted as Xj jG and defined as: 

M = ;U-^(0)/G' 
23 



Its symplectic form um is defined tlirough the relation: 

where p : /U~^(0) — > M is the projection and i : fi~^{0) — > X is the inclusion. 

Let us assume that an integrable Hamiltonian system is defined on X. Let K = 
{Ki, . . . , K^}, X = |dimX denote the set of its integrals of motion. Suppose that this 
system is equivariant with respect to the action of G. This is equivalent to the statement, 
that Ki and /i" form a closed algebra /C with respect to the Poisson brackets. Let us 
assume that on the zero level of the moment map /i the center Z(/C) of the algebra /C is 
sufficiently big, i.e. the dimension of its spectrum equals half the dimension of M. Then 
the integrable system on X descends to the integrable system on M, /C being replaced by 
Z(/C). 

Now let us impose one further restriction. Suppose that X possesses another G- 
equivariant integrable Hamiltonian system, with integrals Q = {Qi,...,^^}, which is 
dual to the system K (algebraically it means that K and Q generate all functions on X). 
We also assume that Q descends to M. 

On the original manifold X the evolution of the system K looks non-trivially in the 
action-angle variables for the system Q and vice versa. The same is true for the reduced 
systems. The advantage of the consideration of X is that the systems on X can be much 
simpler then those on M. In the following sections we shall consider various examples of 
this situation. 

The similar statements hold in the case of Poisson manifolds, the relevant reduction 
being the Poisson one (one first takes a quotient with respect to the group and then picks 
out a symplectic leaf). We leave the details to the interested readers. 

Now we proceed to the explicit constructions. We will discuss the models introduced 
in the previous section on case-by-case basis and show how the reduction which yields 
these systems also explains the dualities between different systems. 

Rational CM model. This model can be obtained as a result of Hamiltonian reduction 
applied to T*g x (9 [^ for g = su{N), O = (DP ~ . The symplectic form on this manifold 



is the sum of Liouville form on T*g and —Nux Fubini-Study form on O. Let (ei : . . . : e^) 
be the homogeneous coordinates on O. The group G = SU{N) acts on T*g via conjugation 
and on O in a standard way (O = G/H, H = S{U{N - 1) x U{1))). Then the moment 
map for the action of ^ on T*g x (9 is 

li = ad*Q{P)-J J,j = iy{N5,j - e^e*) (4.12) 
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where Q G g, P G g*. Now we choose two sets of Hamiltonians: 

Hk = \txP^ and HP = \ttQ^ (4.13) 

k k 

If we identify g* and g with the help of Tr then the equation fj, = has the form: 

[P,Q] = J (4.14) 

which is obviously preserved by the involution: P ^ Q,Q —>■ —P. So we are guaranteed 
to get a self-dual system. Now we have to find suitable coordinates and action variables. 
Let us choose the gauge (remember that we have to mod out ( [4.14| ) by the action of G) : 

Q = diag{qi,...,qN) (4.15) 

This gauge is preserved by the action of the maximal torus T = U{1)^~^ which turns out 



to be sufficient to set all e^ to be equal: e^ = 1 [54]. Then the equation ( 4.14 ) fixes P 
which turns out to be nothing but L in (|4.2|) . As it is obvious that the reduced symplectic 
form equals J2i ^Pi /\dqi (with the constraint J2 ^i — J2 Pi — ^) ^^^ concludes that g^'s are 
the action variables for the system generated by Hj^^s. Therefore eigenvalues of P are the 
action variables for the fiows generated by -ffjt's. We therefore proved the following 
Statement. Consider the map: 

cr■.{{p^,q^)}^m,-V^)} (4-16) 

where rji 's are the eigenvalues of L = P and ^i are the diagonal entries of Q in the 
eigenbasis of P. It is an involution. 

Let us go back to the systems ( |4.13|) . The moment map equation ( [4.14| ) is obviously 



preserved by the transformations of the form 

{P,Q)^{aP + bQ,cP + dQ) ad-bc = l (4.17) 

which form S'L2(IR) group. The transformed Hamiltonians 

g-Hk = ^TT{aP + bQ)'' 



are easy to express through the original Hamiltonians ( [4.13|) in the coordinates [pi, qi) 



g ■ Hk{pi, qi)\„ = Hk{api + bqi, qi 
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Let us restrict our attention to the SL2(2) subgroup of the group ( |4.17|) . It is generated 
by the transformations 

It is clear that S coincides with the involution leading to ( [4.16| ) while T is the unit time 

evolution with respect to the Hamiltonian H^ . 

Trigonometric CM, Rational RS. The trigonometric CM system can be obtained as 



Hamiltonian reduction applied to either T*G x O ||53|| or T*g x O ||55|| where g is the 
central extension of the loop algebra. In the latter case one has to specify the action of 
the gauge group LG on the orbit O. The correct choice is the most natural one: since the 
orbit is finite-dimensional, the only sensible way the loop group can act on it is through 
the evaluation at some point. The elements of T*g of our interest are the pairs: 

P{x),kd^ + Q{x) (4.19) 

where /c is a fixed number, P{x) is a g-valued function on a circle S^ and Q{x) is a gauge 
field on a circle. The phase space is acted on by the gauge group: 

P{x) ^ g{x)-^P{x)g{x), Q{x) ^ g{x)-^Q{x)g{x) + kg-\x)d^g{x) (4.20) 

The moment equation has the form: 

kd^P+[Q,P] = J5{x) (4.21) 



where J is the one from (|4.12|) . The number k can be rescaled by the choice of the radius 
of a circle S"*". Instead we choose the circle of unit radius and keep k. To solve the equation 
( [4.21|) we fix a gauge ( |4.20|) . We can either decide that Q is a constant diagonal matrix 

Q = diag(gi,...,gAr) 

and then the solution for P{x) will produce the Lax operator ( |4.4|) of the Sutherland model 
withi?= ^ il],|5|]. 

It is quite amusing that the same reduction yields the rational RS model as well. In 
order to see that choose the gauge 

P{x) = diag (pi {x),...,pN (x)) (4.22) 
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Then the moment equation ( 4.21 ) imphes that Q{x) is diagonal everywhere except x = 
where it has an off-diagonal part proportional to the delta-function. At the same time 
P{x) is forced to be constant Pi{x) = qi. 

Q{x),j = e,{x)5,j + 5{x)^^ (4.23) 

Qi Qj 

The natural candidate for a Hamiltonian in this setting would be a gauge invariant function 
of Q{x). Since Q{x) is actually a gauge field the gauge invariant function is a trace in some 
representation R of a Wilson loop: 

Hr = Tr^Pexp / ^Q{x)dx (4.24) 

which is easy to evaluate provided we assume the following structure of the diagonal piece 
of Q{x) (which supported by the alternative derivation of the solution to the moment 
equation below): 

9i{x) = ipi{x) +5{x)^ 



lU 



fc^. ^^ - ^» 



which makes the Wilson loop 



S = Pexp j^Q{x)dx = diag (^e^ /^»^"^''") exp f'-^r\ (4.25) 



with r being the matrix: 



^ij 



— , i^j, rii = - J^ Yij (4.26) 

Qij 



i/i 



It is shown in the Appendix (in the trigonometric case from which this one follows as 
well) that the matrix B is gauge equivalent to (|4.8| ) with the identification (3 = ^, and (cf. 
(H)): 

Pi + '^^^^^^^t/(P7) = f ^tix)dx. 

One can also get the same matrix (without the assumptions like ([4.23|) ) B by perform- 
ing a reduction of T*G under the adjoint action of G at the same level J of the moment 
map: 

l^ = B-^PB-P = J 

Trigonometric RS = Relativistic trigonometric CM model. There are three different 
approaches, all leading to the same finite-dimensional Hamitlonian system. There are two 
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Hamiltonian reductions and one Poisson reduction. The advantage of Hamiltonian one is 
the simphcity and geometric clarity. The advantage of Poisson one is finite-dimensionahty 
at each step and considerable simplicity of the proof of the canonical commutation relations. 
We try to outline all three approaches with the emphasis on the Poisson reduction, as the 



relevant Hamiltonian reduction was described in some details in ||5^. We keep in mind a 
sequence of contractions: 

At^ -^ T*G ^GxG (4.27) 

where the first entry is the space of G- valued gauge fields on a two-torus T^, the second 
entry is the cotangent bundle to the central extension of the loop group LG and the last 
one is the space of lattice (for the simplest graph, representing a two-torus) connections, 
described below. 

Hamiltonian approach. Consider the space At^ of SU{N) gauge fields A on a two-torus 
T^ = S^ X S^. Let the circumferences of the circles be R and (3. The space At^ is acted 
on by a gauge group Q , which preserves a symplectic form 

O = -^ /" Tr5A A 5 A, (4.28) 

47r'^ J 

with k being an arbitrary real number for now. The gauge group acts via evaluation at some 
point p e T^ on any coadjoint orbit O of G, in particular, on O = CP'^"^. Let CP'^"^ 
have a —Nux Fubini-Study symplectic form. Let (ei : ... : cat) be the homogeneous 
coordinates on O. Then the moment map for the action of Q on At^2 x O is 

kFA + J6\p), J,j = w{5,j - e,e*) (4.29) 

Fa being the curvature two-form. Here we think of e^ as being the coordinates on C 
constrained so that Yli kiP = ^ ^^^1 considered up to the multiplication by a common 
phase factor. 

Let us provide a certain amount of commuting Hamiltonians. Obviously, the eigen- 
values of the monodromy of A along any fixed loop on T^ commute with themselves. We 
consider the reduction at the zero level of the moment map. We have at least A^ — 1 
functionally independent commuting functions on the reduced phase space M.^- 

Let us estimate the dimension oi M.jy. If z^ = then the moment equation forces the 
connection to be flat and therefore its gauge orbits are parameterized by the conjugacy 
classes of the monodromies around two non-contractible cycles on T^: A and B. Since 
the fundamental group 7ri(T^) of T^ is abelian A and B are to commute. Hence they are 

28 



simultaneously diagonalizable, which makes A^o a 2{N -- 1) dimensional manifold. Notice 
that the generic point on the quotient space has a non-trivial stabilizer, isomorphic to the 
maximal torus T of SU{N). Now, in the presence of O the moment equation implies that 
the connection A is flat outside of p and has a non-trivial monodromy around p. Thus: 

ABA-^B-^ = exp(i?/3J) (4.30) 

(the factor i?/3 comes from the normalization of the delta- function in ( ^.29| )). If we diago- 
nalize A, then B is uniquely reconstructed up to the right multiplication by the elements 
of T. The potential degrees of freedom in J are "eaten" up by the former stabilizer T of a 
flat connection: if we conjugate both A and B by an element t & T then J gets conjugated. 
Now, it is important that O has dimension 2(A^ — 1). The reduction of O with respect to T 
consists of a point and does not contribute to the dimension of A^ j, . Thereby we expect to 
get an integrable system. Without doing any computations we already know that we get a 
pair of dual systems. Indeed, we may choose as the set of coordinates the eigen- values of A 
or the eigen- values of B. The logarithms of the eigen- values of B are the action variables 
for the system generated by TrB^, and vice versa. 

The two-dimensional picture has the advantage that the geometry of the problem 
suggests the S'L2(2)-like duality. Consider the operations 5" and T realized as: 

S:{A,B)^{ABA-\A-^); T : {A, B) ^ {A, BA) (4.31) 

which correspond to the freedom of choice of generators in the fundamental group of a two- 
torus. Notice that both S and T preserve the commutator ABA~^B~^ and commute with 
the action of the gauge group. The group F generated by S and T (it is a subrogup of the 
group OutPree(2) of the outer authomorphismes of the free group with two generators) 
seems to be larger then SL2(2). However in the limit P,R —>■ it contracts to SL2(2) in 
a sense that we get the transformations ( [4.18|) by expanding 

A = l + pP + ..., B = 1 + RQ + ... 

for R,(3^0. 

The disadvantage of the two-dimensional picture us the necessity to keep too many 
redundant degrees of freedom. The flrst of the contractions ( [4.27| ) actually allows to replace 
the space of two dimensional gauge flelds by the cotangent space to the (central extension 
of) loop group: 

T*G = {{g{x),kd^ + P{x))} 
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which is a "deformation" of the phase space of the previous example {Q{x) got promoted 
to a group- valued field). The relation to the two dimensional construction is the following. 
Choose a non-contractible circle S"*^ on the two-torus which does not pass through the 
marked point p. Let x,y he the coordinates on the torus and y = is the equation of the 
S-*^. The periodicity of a; is /3 and that of y is R. Then 

P{x) = A^{x, 0), g{x) = Pexp / Ay{x, y)dy. 

Jo 

The gauge transformations on S-*- transform on {g(x), P{x)) is a way, similar to ( [4.20|) . 
The moment map equation ( [4.29|) goes over to the moment map equation [^ : 



kg-^d^g + g-^Pg -P= J5{x), (4.32) 

with k = ^. The solution of this equation in the gauge P = diag(q'i, . . . , qat) leads to the 
Lax operator A = g{0) of the form ( |4.11| ) with R, P exchanged [^. On the other hand, if 
we follow ( f4.22|) and diagonalize g{x): 

g{x) = diag {zi = e*^«^ ,...,zn = e^^"^^) (4.33) 

then a similar calculation leads to the Lax operator 

B = Pexp (h -P{x)dx = diag(e'^*) expzi?/3i 



WY 



with 

thereby establishing the duality A ^^ B explicitly. 

Poisson description. Here we introduce a set of commuting functions on the space of 
graph connection on a graph, corresponding to a moduli space of fiat connections on a 
torus with one hole and describe the flow generated by this set. Being reduced to a 
particular symplectic leaf of the moduli space of flat connections on the torus , this set 
of functions turns out to be a full set of commuting Hamiltonians. We introduce another 
full set of commuting variables and write down the Hamiltonians taking the latter set as 
a set of coordinates thus recovering the Ruijsenaars integrable system. Consider a graph, 
consisting of two edges and one vertex with the fat graph structure corresponding to a 



punctured torus [^ . The space of graph connections A for such graph is just a product 
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of two copies of the group G: A^ = G x G = {{A, B)\A, B G G}, where A and B are 
assigned to the edges of the graph. For a choice of cihation on A^ the Poisson bracket on 



A is given by the relations, foUowing from the general rules [57 



{A^A} = r"A(g)A + A® Ar"-2(A® l)r"(l® A) 

{B<^B} = r"S®S + S®Sr"-2(S®l)r"(l(g)S) (4.34) 

{A^B} = r{A(^B)+A(^Br+{l0B)r2i{A®l)-{A0l)r{l0B), 

where r" = ^{r — r2i). 

Now let us restrict ourselves to the case G = SLn and the standard r-matrix: 

r = ^ £;, ® £;_, + ^ ^ if, ® if,, r" = ^ ^ £;, A e_^ (4.35) 

a>0 i a>0 

In this case one can easily derive the following commutaion relations 

{TrA'^, A} = {TrS'", B} = (4.36) 

{TrA", B} = n(A")o {TrS", A} = nA(S")o (4.37) 

where (^)o denotes the traceless part of the matrix X . Therefore, the functions Tri?" for 
n = 1 . . . A^ — 1 considered as Hamiltonians generate commuting flows on A^ . 

S(ti,...,t^_i)=S(0,...,0) 

A (ti, . . . , tA,_i) = A (0, . . . , 0) el*^^+-+*--^ )o 



As it was shown in [^ the lattice gauge group Q^ acts on Q^ in a Poisson way, and the 
quotient Poisson manifold coincides with the moduli space A4 of smooth flat connection on 
the Riemann surface, corresponding to the fat graph L. In our case the group Q^ is G itself 
(for the graph has just one vertex) which acts on A and B by simultaneous conjugation. 

g:{A,B)^{gAg-\gBg-'). (4.39) 

The functions TrA^ and Tri?'^ are invariant under this action, and therefore their 
pull-downs on the moduli space Ai generate commuting flows there, which trajectories are 
just projections of (|4.38| ). 

However the moduli space Ai in our case is a Poisson manifold with degenerate Poisson 
bracket. The Casimir functions of this Poisson structure are the functions of conjugacy 
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classes of monodromies around holes and constant value levels of such functions are just 
the symplectic leaves of M.. In our case such Casimir functions are Tt(^ABA~^B~^^ , 
pulled down to M.. 

Different symplectic leaves have different dimensions and the lowest dimension of them 
is 2(A^ — 1). These leaves correspond to the monodromy around the hole conjugated to a 
matrix 

rkP < 1, z^ is a numerical constant from the previous section parameterizing the set of 
symplectic leaves of lowest dimension. Let t = e~'^^^'^ . On the leaf TWjy the family of 
functions TrA'^, k = 1, . . . , N — 1 forms a full set of Poisson-commuting variables. 

Introduce local coordinates on these symplectic leaves in the following way. Let 
zi = e^^'^^ , . . . , zn = e^^'^^ be the eigenvalues of the operator A and ^i, . . . , fiN are 
the corresponding diagonal matrix elements of B (in the basis, diagonalizing A). One can 
check that in this basis 

Si = VJ^jj^^- (4.40) 

The functions Zi and fij are well-defined locally on the symplectic leaf A4^. Their Poisson 
brackets are equal to: 

{zi.Zj} = 



X^ij f^j f ZifXjOi^ 



J- 



To define the variables, canonically conjugated to Zi we can just multiply fXi by factors 
independent on /i^. For example one can take: 



-'--'- \/ [zk -tZi){Zi -tZk) 



One can check, that these new variables s,- have the Poisson brackets 



^i 



{si, Sj} = {zi, Sj} = ZiSjSij. (4.43) 

Substituting this back to the formula (|4.40|) we get: 



S., = -i-^ ($+$-$+$-) ^/^ (4.44) 
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which is gauge equivalent to (|4.11j ). 

Moral revisited. We have seen in all the previous examples that the origin of the dual 
system is connected with the existence of transversal G-invariant foliations on the original 
space, which become Lagrangian foliations when pulled down to the quotient. The sim- 
plicity of the operating with dual systems in the advocated framework in the classical case 
allows one to hope that the duality can systematically elevated to the quantum case as 
weU. See [|8l||. 



4.3. Appendix. Computation of the Poisson brackets 

The bivector defining the Poisson structure on A^ can be rewritten in the form 

^= ^ E Ep®El^-\e{u,v)+e{i,3)), (4.45) 

where e(z, j) is —1, or 1 depending on whether i is less, equal or greater than j respectively 
and E^ are the standard GLtv generators acting on the u-th. end of the edge. (In our 
case i?^ acts on A from the left, -E^ acts on B from the left, -E^ acts on A from the 
right and E^ acts on B from the right.) 

It is not convenient to compute the Poisson brackets between Zi and /Xj using this 
bivector directly, for it does not preserve the diagonality of the matrix A. However, to 
compute the Poisson brackets of the gauge invariant functions we are allowed to add to 
this bivector any term vanishing on such functions i.e. any terms of the type ^^ E^ ® X 
or X 'Yliu-Ej \ where X is an arbitrary vector field. Using this and also the fact that 
for the diagonal A: 






Zj 



J 'y . -y . J J 



(4.46) 



z 



J 



the bivector n can be transformed to the form: 



„' = j; £f ) A £* • ^i±^ + i 5; Sf A £■'" + Ef> A £f > (4.47) 

Applying this bivector for the chosen A and B we get the desired Poisson brackets, b 



The equations (|4.46|) are the infinitesimal forms of the statements that (i) up to the 
gauge transformation the conjugation of A hj g is equivalent to the conjugation of B by 
g-^ thanks to (|]3|) ; 
(ii) gLA = A-\Ag^^)A, with gL = g^^ 
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4-4- Appendix. Solution of the moment equation 
Here we solve the equation ( ^4.301) : 



A-^BAB-^ =expi?/3J 

(4.48) 
J = -ii/{Id-e(g)e^), {e\e)=N 



with A, B - N X N unitary matrices defined up to the gauge transformations (|4.39| ). We 
use the notation: a = R(3iy. We partiaUy fix a gauge: 

A = diag {e'^'J' , . . . , e^^'^^ ) (4.49) 

which leaves gauge transformations of the form 

h = exp {i diag(/i, . . . , In)) . (4.50) 

which preserve A, conjugate B and map e to h~^e. The exponent expi?/3J is easy to 
compute: 

exp J = e"'" ( Id H e (g) e^ 

Let / = B-^e, z^ = e'^'i\ $+ := |e,|2, $" := |/,|2. Then: 

ANa 






(4.51) 






i=l ^ 

The last equation implies (see below): 



^— Z_/^'t PW = n('-..) (4.52) 



i=l 



Now the unitarity of B implies, that 



^^'^ = ^^^^* ^^ ^5 ^.i^'(-.)(^.A.-e-)(^.A,-e— ) ^^-^^^ 

Hence 

N P(e''°'z) 
$- = - ^^ ^^^ (4 54) 



"I-' V^lyl 
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To prove ( 4.52 ) consider the contour integral 

1 /• P{e-''^z)dz 



2ni /]r P{z){z- e'-'^Zj] 
To prove (|4.54| ) consider the integral: 



ReSr 



-iNa 



1 f P{e-''^z)dz _ _ P'jzk) 

2m fr P{z){z - e-z.)(^fc - e— ^) '^^^^'^'"^'= "'' P{e^-z,) 

In both cases the contour JT surrounds the roots of P{z). 
Notice that both $„ are real: 



^^ _ iVsin(«/2) -p. sin 



Rqij±a 



sm{Na/2)j^^ sin^^ 



(4.55) 



Substituting this back to ( ^4.51| ) we get: 



B 



1-3 



^i{{l-N)a/2+Rqi, /2+Si-ipj) _ 



-m 



sm 



A^'sin 



Rqij+a 



$+$- 



(4.56) 



where e^ =: |ei|e*^% /j =: |/j|e*'^^. The gauge transformations ( ^4.501 ) allow us to set 
(/?, + Rq^/2 = 0. Then define 



P^ = -^ ((1 - N)a/2 + Rqj2 + £,) 
Finally, the matrix B can also be written as: 

B = ($")"^ {e-'^Pe-'""') ($-)^ 
where $~ = diag($~), p = diag(pi), 



a + n i / $„• 

Pi=Pi 777^ + 777^m 



(4.57) 



(4.58) 



2(3 2(3 ^V$+ 



(4.59) 



and 



Zi 



<-ij 



Zi Zj 



-, i^j 



1 z,P"{z,) 

2 P'{zi) 
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(4.60) 



To prove the last statement consider the matrix 



Rij{a) 



- (^) .+ 



sm 



i(N-l)a Z 



A^sin 



Rqij+a 




Pie-'^^Zi) 



z^ (e-^^z, - Zj)P'{zi) 



We have: 

B = ($")"^ {e-"^PR{a)) ($-)^ 

A simple contour integral calculation shows that 

R{ai)R{a2) = R{ai + 02) 

The rest follows by expanding near a = 0. If one performs an expansion near R = one 
gets the statement that the rational Lax operator ( ^7^ ) is conjugated to the operator of 
the form announced in (|4.25| ). 

It is amusing that the expression 2iglog($^/$~) appears quite often in Bethe Ansatz 



Equations for XXX magnets and their field theoretic limits [59 . 
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5. Gauge theories and duality in integrable systems 

5.1. Old approach: many-body systems as low- dimensional gauge theories 

It is a fruitful approach to think of the many-body system as of the gauge theory of 
a certain kind. Namely, the particles of the model can be identified (sometimes) with the 
eigenvalues of the Wilson loops in the theory and the gauge dynamics becomes a dynamics 
of the particles. Of course, in the real four dimensional world the gauge field has infinitely 
many degrees of freedom and we don't expect to see any tractable quantum mechanical 
system unless we have a principle which allows us to restrict the dynamical problem to 
a finite number of degrees of freedom. The simplest case is the case of low-dimensional 
gauge theory, where the gauge field simply doesn't have propagating degrees of freedom. 
Consider, for example, two dimensional Yang- Mills theory with a gauge group G = U{N). 
When formulated on a circle of radius R in the Hamiltonian formalism the theory has as 
a phase space the space of gauge fields A{x) on the circle and their duals - chromoelectric 
fields E{x). The gauge group acts on {E, A) as follows: 

(E, A) ^ {g-'Eg, g-'d^g + g-'Ag) (5.1) 

leading to the Gauss law dxE -|- [A, E], which is nothing but the moment map from the 
section 4.2. We can go to the gauge where A is a constant (w.r.t. x) diagonal matrix 

A = diag{qi,...,qN) (5.2) 

Here are our particles. The time evolution makes qi to move and depending on the cir- 
cumstances such as the presence of the sources like J (which correspond to the time-like 
Wilson lines) one gets the Hamiltionian system of the kind we described and studied. 
The large gauge transformations shift g^'s by integer multiples of i- making them live 
on a circle of radius 2^- ^^^ '^^^ S^^ more complicated examples by deforming the 
model as follows. Replace S-*^ by T"^, E by the second component of the gauge field 
along the torus, the symplectic form being L,2 Tr5A A 6A. Then the Gauss law becomes 
Fa = dxAy — dyAx + [A^, Ay]. Setting it to zero allows to diagonalize A^, Ay simultane- 
ously: 

(t)^-<(::)--fe)) '-' 

Here, Xi and yi do not Poisson-commute, although both live on circle. One gets, therefore 
a system of relativistic partcles on a circle. The radius of the circle is ^J^ , the speed of 
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light is Ry. The gauge theory this model corresponds to is known as Chern-Simons theory 
on a torus (perhaps with punctures). 

One can go higher in dimensions with some care. For example, by considering super- 
symmetric M = 2 theory in (i = 4 with compact space M"^ one gets a quantum mechanics 
on the moduli space of monopoles in IR . 

5.2. New approach: many-body systems in supersymmetric gauge theories 

Recent progress in the understanding of non-perturbative phenomena emerged after 
the work of Seiberg and Witten on four dimensional M = 2 SYM [jl2| and works of Seiberg 
and his collaborators on A/" = 1 (i = 4 theories. The major tool in these studies is the 
low energy effective Lagrangian which is constrained by two priniciples - the holomorphy 
of chiral objects and electric-magnetic duality. It is the electric-magnetic duality which 
makes the integrable systems to appear in the solutions to the gauge theories. 

In particular, one can argue on the general grounds [0 that any A/" = 2 supersym- 
metric gauge theory in four dimensions corresponds to a certain integrable system in the 
holomorphic sense. The point is that the Coulomb branch of the theory parameterizes 
the family of abelian varieties (whose period matrix coincides with the matrix of coupling 
constants of the effective low-energy abelian theory) . Moreover the total space must carry 
a holomorphic symplectic form uj, whose integral along the cycle in the fiber gives rise to 
a derivative of the central charge of a BPS representation of A/" = 2 susy algebra along 
the base. Moreover the abelian varieties must be Lagrangian with respect to uj. 

The integrable systems corresponding to a large number of field theories are identified. 
In particular, the low-energy theory of the pure N = 2 SU{Nc) SYM is governed by 
the An^-i periodic (or affine) Toda system. The M = 2 theory with a massive adjoint 
hypermultiplet corresponds to the elliptic Calogero-Moser system, where the mass (which 
is naturally a complex parameter in the M = 2 theory) is identified with the coupling 
constant. The theory is UV finite (in fact, it is softly broken A/" = 4 theory) and therefore 
has as another modulus - the ultra-violet coupling r which enters the integrable model as 
the modulus of the curve. Another theories which were mentioned so far are the relativistic 
generalizations of those two. These correspond to five dimensional gauge theories with the 
same number of supercharges, compactified on a circle of a finite radius R. The speed of 
light of the relativistic model is proportional to the inverse radius -^ of the circle. For the 
theories with fundamental matter the firm identification with the integrable systems has 



been made in four |j2^ [^^ as well as in five and six dimensions [^ 
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In some cases the dualities suggested by the integrable systems are not obvious on 
the field theory side. We plan to return to more detailed treatment of these cases (which 
involve six dimensional theories) in the future. 

5.3. Dualities in field theories vs. dualities in many-body systems 

Dualities in the old approach. Let us start with the two-dimensional Yang-Mills the- 
ory with the gauge coupling g'^ formulated on a Riemann surface of area A. It was shown 
by E. Witten in pTI] that the perturbative in g'^A part of the correlation functions in this 



(non-supersymmetric) theory coincides with the correlation functions of certain observables 
in twisted A/" = 2 supersymmetric two-dimensional Yang- Mills theory. 

Among the twisted supercharges of the latter theory one finds a scalar Q which annihi- 
lates the complex scalar (p in the vector multiplet. The observables constructed out of the 
gauge invariant functions of (j) and their descendants can be mapped to certain observables 
in non-supersymmetric theory. 

As we discussed above, when Yang-Mills theory is formulated on a cylinder with the 
insertion of an appropriate time-like Wilson line, it is equivalent to the Sutherland model 
describing a collection of A^ particles on a circle. The observables Tr^'^ of the previous 
paragraph are precisely the integrals of motion of this system. 

One can look at other supercharges as well. In particular, when the theory is for- 
mulated on a cylinder there is another class of observables annihilated by a supercharge. 
One can arrange the combination of supercharges which will annihilate the Wilson loop 
operator. By repeating the procedure similar to the one in |^T[] one arrives at the quantum 
mechanical theory whose Hamiltonians are generated by the spatial Wilson loops. This 
model is nothing but the rational Ruijsenaars- Schneider many-body system. 

The duality between these two systems is a consequence of the fact that when lifted 
to the supersymmetric model both field theories become equivalent to the same J\f — 2 
super- Yang-Mills theory in two dimensions. 

The self-duality of trigonometric Ruijsenaars system has even more transparent physi- 
cal meaning. Namely, the field theory whose quantum mechanical avatar is the Ruijsenaars 
system is three dimensional Chern-Simons theory on T^ x R^ with the insertion of an ap- 
propriate temporal Wilson line and spatial Wilson loop. It is the freedom to place the 
latter which leads to several equivalent theories. The group of (self-) dualities of this model 
is very big and is generated by the transformations S and T ( |4.31j ). 
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In short, the duahty reveals here itself as a consequence of Lorentz invariance of the 
underlying field theory. 

Duality in the new approach. The new approach deals with supersymmetric gauge 
theories in three, four, five and six dimensions. Perhaps the richest case is the six dimen- 
sional theory compactified on a three dimensional torus T"^ down to three dimensions. 

As was discussed extensively in [jl5| in case where two out of three radii of T"^ are much 
smaller then the third one R the effective three dimensional theory is a sigma model with 
the target space X being the hyper-kahler manifold (in particular, holomorphic symplectic) 
which is a total space of algebraic integrable system. The complex structure in which X is 
the algebraic integrable system is independent of the radius R while the Kahler structure 
depends on R in such a way that the Kahler class of the abelian fiber is proportional to 
1/R. 

The duality of the integrable systems shows up in the gauge theories in the several 
ways. 

First of all AA duality the well-known phenomenon in the four dimensional M = 2 



gauge theory which was observed and exploited in [|T2[ and then later on in the plenty of 
works. The low-energy effective theory has different sets of relevant degrees of freedom 
over different regions of the moduli space of vacua. The transformations between differ- 
ent descriptions go through the electric-magnetic duality on the gauge field side which is 
accompanied by supersymmetry by a AA-type duality on the scalar side. Although this 
duality is connected with the electric-magnetic symmetry which is not realized geometri- 
cally in four dimensions, it does become geometric when the theory is lifted to a tensor 



theory in six dimensions |62 



The duality of the AC type is also present and is rather interesting. As one varies the 
moduli of T"^ the geometry of X varies as well. In particular, different four dimensional 
theories can flow to the same three dimensional theory. This is where the AC duality in 
the integrable systems shows up. 

For example, a certain scaling limit of the five dimensional SU{N) theory with massive 
adjoint hypermultiplet, compactified on a circle seems to be equivalent /dual to four di- 
mensional SU{N) theory with massive adjoint hypermultiplet when instanton corrections 
are turned off in both theories. 

The theory in three dimensions which came from four dimensions upon a compact- 
ification on a circle whose low-energy effective action describes only abelian degrees of 
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freedom can be always dualized to the theory of scalars/spinors only, due to the vector- 
scalar duality in three dimensions. In this way different sets of vector and hypermultiplets 
in four dimensions can lead to the same three dimensional theory (one of the examples of 
such symmetries is provided by the three dimensional mirror symmetry |^^ [Q ) . 



One can also use the AC duality to establish the following fact. Consider two AC dual 
integrable systems. Consider three dimensional A/" = 4 supersymmetric gauge theory whose 
Higgs branch is Af - the phase space of these systems. Then two M = 2 supersymmetric 
theories whose spaces of scalars are both X but the superpotentials are taken from the 
sets of Hamiltonians of the first and the second systems respectively are dual to each other 
in the sense that both flow in the UV/IR (depending on the whether these Hamiltonians 
correspond to the relevant or irrelevant operators) to the same theory. 

We are certain that there are more applications of the notion of duality in integrable 
systems both in the theory of integrability itself and in the physics, gauge theories being 
the arena for the most immediate ones. 
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